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ABSTRACT: In this paper we present the revised view on the optical activity 
of crystals based on the model of two dumped coupled oscillators. The results are 
compared with the results of the same problem solved before and it is presented 
that the results of the quantum mechanical model are the generalization of the 
classical model only and they do not introduce the new terms of quantum me- 
chanical nature. The results are discussed with regard to the rotational strengths 
of normal modes of vibrations and the formula for the complex rotatory power 
p{ui) containing the rotational strengths is presented. 

1 Introduction 

Optical activity (OA) is manifested by dissymmetric species because the tran- 
sition probabilities are different for the left and the right circularly polarized wave. 
This means a difference between the complex refractive indices that we denote fii 
for the left and n r for the right circularly polarized wave. The complex rotatory 
power may be considered to be complex quantity given by 

p = p + ia = — {n x - n r ) = — [(ni -n r )+i (kj - n r )} , (1) 

where p is the rotation angle of incidenting linear polarized light per unit length, 
and for a very excellent approximation a is the ellipticity per unit length. The 
variation of p with frequency uj or wavelength A is called the optical rotatory 
dispersion (ORD) and the variation of a is the circular dichroism (CD); ni, n r 
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and tv[ , Kjj> £1X6 real and imaginary parts of the complex refraction indices ni 
and n r . So that the OA has both aspects arising from interaction of radiation 
with matter - dispersive and absorptive and these aspects are connected by the 
Kramers - Kronig relations. 

The important group of the optically active crystals are the crystals with 
screw axis of symmetry belonging to the space group of symmetry D\ and its 
enantiomorphic D\. The typical representatives of these crystals are a-quartz, 
cinnabar, tellurium, selen, camphor and benzil. The optical activity of these crys- 
tals is caused by the assymetrical originating of the crystal structure because the 
molecules or atoms forming the crystals are symmetrical. Camphor is the excep- 
tion because its molecules are optically active. It means that the other crystals 
belonging to these groups of symmetry are optically active in the crystalline state 
only. 

The dispersion of the optical activity was studied in works based on the ex- 
citons theory |l|, |, ||, |], on the theory of coupled oscillators || ||, |7|, |8|, || or 
on the Lagrangian formalism [TO. But all these models based on the exciton 



theory, some models of coupled oscillators and the model in [Kj are solved in 
the frequency region far from the absorption range. It is well known that the 
CD is nozeroth only in a very narrow frequency region in the absorption range 
and therefore these models solve only ORD as one part of the optical activity 
independently on the CD. 

Both aspects of optical activity can be solved by means of the model of the 
dumped coupled oscillators. Vysin |8j has solved this model semiclassically and 
Janku || by quantum mechanical way where the role of dumping plays the limited 
lifetime of oscillators in excited states. But Janku has got some other terms 
in comparison with the results of || in his formulae for ORD and CD and he 
supposes that these terms are of quantum mechanical nature. At the same time 
the sense of these terms is not evident. For this reason we revise the solution of 
this problem. 

All above mentioned works that solve the optical activity of crystals by the 
coupled oscillator model arised from the Kuhn model |TT| of two coupled oscilla- 
tors forming the compound oscillator. It is well known that two asymmetrically 
oriented harmonic oscillators coupled together form the fundamental optically 
active unit. With regard to the structure of crystals belonging to the space 
groups of symmetry D\ and D\ we can use the following Chandrasekhar model 
of compound oscillator. 

We assume that the coordinate axis z is parallel with the crystal axis c. The 
first single dumped linear harmonic oscillator lies in the plane z = and his 
direction of vibrations is given by direction cosines a, /3, 7. The second oscillator 
lies in the plane z = d where as d we denote the projection of distance between 
oscillators into the z axis and his direction of vibrations is turned by the angle 
9 around axis z with respect to the first one. The direction of vibrations of the 
second oscillator is so given by the cosines a cos 9 — j3 sin 9, a sin 9 + (3 cos 9, 7. 
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Both oscillators lie on the helix which is given by the structure of crystal. The 
angle 9 is, of course, 120 deg for the crystals belonging to the space groups of 
symmetry D\ and D®. 

Both oscillators forming the compound oscillator are identical. They have 
the identical mass, electric charge and lifetime in excited states. We take the 
interaction between adjacent oscillators as week, that is for example as dipol - 
dipol character. We suppose that the left and the right circularly polarized wave, 
into which the incidenting linear polarized wave is decomposed in the optically 
active crystal, propagate only along the crystal axis. This is the most practically 
important case because the optical activity is not masked by birefringence. 



2 The quantum mechanical solution of the model 

We can set the Schrodinger equation for the one compound oscillator in the 
field of the left and the right circularly polarized wave 

** = ± |f + ^ ± ^ + ^ + JL ± (2) 

where r±, r 2 are the displacements of oscillators from equilibrium, \1r\r2 is the 
potential energy of mutual interactions of the oscillators, F l { r = eE 1 / ', F l 2 ' r = eE l 2 ' r 
are the electric forces projections of the left and the right circularly polarized wave 
into the vibration directions of the oscillators, e is the charge of electron. The 
upper index l,r holds for the left and the right circularly polarized wave, pi and 
P2 are the moment operators of both oscillators. The small positive parameter 
7o gives the possibility of the adiabatic interaction at the time t = — 00. The 
dumping due to the limited lifetime of oscillators in their excited states is formally 
introduced by the parameter 70. 

The electric vector E of the left and the right circularly polarized wave prop- 
agating along z axis has the components 

E l x = E e- i(wt ~ k ' z) , E\ = Eoe-'i"*-* 1 *-*) 

E r x = E^ l{wt ~ Kz \ E r y = E e"< wt - krZ+ ^ (3) 

and the vector E projections of these waves into the directions of vibrations of 
oscillators are 



E l { r = Eoe-™ (a + 



E l i r = E e~ iuJt 



(a cos9-P sin 9) + (a sin 9 + (3 cos 9) e'K^*) 



, (4) 



where <pi >r = k^ r = n^ r ud/c is a complex phase shift and n^ r are refractive indices 
of medium for the left (index /) and the right (index r) circularly polarized wave. 
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In all terms in eq. (^) and further the + sign holds for the left and the — sign 
for the right circularly polarized wave. 

Due to the mutual interaction between coupled oscillators the natural fre- 
quency ujq splits into two adjacent frequencies of the normal modes of vibrations. 
Introducing the normal coordinates q\ and q 2 which are given by the relations 



n 



1 

V2 



(Qi + Q2) 



(5) 



into eq. (g) we get 
WH = 

+ 



h 2 2 <9 2 ^ moo 2 2 2 _ 1/2 2\ t 



77=1 



±(Fy + F^)p qi + ±(Ft r -F>>' 



Pq2 



e 70 '*. (6) 



It may be easily verified that in the eq. (f|) the expresions (fI' t + F 2 
and -j= {F l { r — F 2 j are the projections of electric forces in normal coordinates 



We denote them as F q f and F*£. 



Now the equation (|6]) can be separated by means of 

^x(qx,t)y 2 {q 2 ,t), 

Tl.q 1 + 7~lq 2 . 



n 



(7) 



Substituting (0) into (H), further dividing by ^i(qi, t)^/ 2 (q 2 , t) and posing the 
member with q 1 to H qi and q 2 to H q2 we get two equations 



h2 ^^K ^^ glM^t) 



2m 

i 



-,2 a2 



dq\ 2 
K%)e Toi fi(9i,i) 



h 2 d 2 ty 2 (q 2 ,t) mu 2 - // 



2m 

i 

muj 



+ 



dql 2 
(F l q2 r p q2 )e^ 2 (q 2 ,t) 



: g 2 2 ^2(g 2 ,t) 



(8) 



and we see that the natural frequency is split into two vibrations of the normal 
modes. For these frequencies we have 



ujf = to 2 + Q, 



UJn 



ujI - Q; 



(9) 



Q = fi/m. 

For F l q ^ and F^ in the equations @ we can derive using eq. (£|) the expression 



77 = 1,2. 



(10) 



The coefficients a l £ are determined by the relations 



(a£) = (a 2 + p 2 ) (1 + cos 0=F0/, r sin0), 

«) 2 = ( a2 + / ?2 ) (l-cos0±0,, r sin0); (11) 



cr^ r are the meaning of the phase shifts only. In the eqs. ( p]) the upper sign 
holds for the left and the lower sign for the right circularly polarized wave. The 
eqs. (0) can be rewritten in the general form 

Hq^n^t) - -— — 2 + -^-q^ v {q v ,t) 

+£; «) E p qv e-<" t+ ^) + ^%(q v , t); (12) 
7] — 1,2. The Schrodinger equations for the normal modes of vibrations then are 

ih d%^ = i^y vM + #£ %M (13) 

where is a nonperturbed and a perturbed hamiltonian. In our case 



and we see that the nonperturbed hamiltonian depends only on the index 77 of 
the normal mode of vibrations and it doesn't depend on the polarization of the 
light wave. We can further write Ti 1 ^ = 7i.g . On the other hand the perturbed 
hamiltonian 

wr-^K;)^^*)^' (15) 

depends on the index of the mode and also on the polarization of the wave. 

The mean value of the induced electric dipol moment from the side of the left 
and the right circularly polarized wave that we hold as small perturbation we can 



solve by means of the Kubo theorem [FJ 

,2 



(*) *) = <« («) *>„ + £ «)* * {Ml e-'i^O^'. (16) 



The first term on the right side of eq. (P!|) is the constant dipole moment 
of system. This term has no meaning for us. The expression \\qn)Pq v *)/ is the 

Fourier transform of the retarded Green function \(q v ,Pq v }/ of operators q v and 
p qv , that is 

00 

W*,» tf =ft / e^ ot ((q v ,p qv )) t dt. (17) 
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For the retarded Green function \(QmPq v )/ m the ground states of quantum 
system \r]o) we hold the relation 

t = -^(DiVol %{t),Pq v ] \Vo)- (18) 

In the eq. fll8D is the unit step function, q n (t) is the operator of q v in the 
interactions representation 

l v ( t ) = e^'Xe-^* (19) 
and for the aplication of the operator in the exponent on the wave function holds 

e ±m^t/h^ = e ±iE nn t/n^ = e ±^n^ n ). ( 2 0) 

In following rearangement we use that for matrix elements of the production 
of operators F q and K q holds in the normal modes the relation 



(Vm I Fq n Kq v \Vn) = (Vm I F q v \Vk) (Vk I Kq n \ T] n ) 
k 

and further we use that the matrix elements of operator p q are 

(Vm\Pq v \Vn) = imuJ Vmn (r) 



(21) 



(22) 



'/mn 'lm '/n 



Then we can solve that 

(Vol %(t),Pq«} \Vo) = tmJ2^ v J(Vk\q v \m)\ 2 ■ (e^°* + e^ *) (23) 



and after substituting the result of eq. ( p3|) into ([18]) and solving the integral on 
the right side of eq. (|17D we have 



W 9 „ 



2ium^ ^nkoKVklqnlVo)? 



Ti 



y 



(24) 



The induced dipole moments d l q ,r we can solve by means of (|24"D and (|16|) 

2e 2 {a l q (f uJ n J{r] k \q n \^)\ 2 



J2 



(25) 



and we can introduce the oscillator strengths of the normal modes of vibrations 
into (E5|) by the relation 



2mw rt J(ri k \q rj \ri )\< 

h 



(26) 
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Now the induced dipole moments can be expressed as 

k "" ~~VkO 

The mean polarizability per unit volume would be then 



< = E«) : 



(27) 



X l ' r 



Nd I,, 



£ oe -i(Wa^)+7oi' 



(28) 



where N is the number of coupled oscillators in a volume unit. But if we include 
all coupling in the crystal then we see that the number of coupling between two 
adjacent oscilators in the direction of propagating light is the same as the number 
of single oscillators. 

Substituting from fl2"g|) and ( |2"T| ) into the Drude - Sellmaier dispersion relation 
of refractive indices of the crystals 



l,r 



v =l 



we have 



^2 



n 



4ttA^ ^ 



and we are able to solve the relation 

„ AirNe 2 2 



n x -n r 



m 



EE 

r)=l k 



«) 2 -K) : 



In 



'/AO 



uS 2 — 2i"f uj 



that after substituting from ( |TT|) gives 
47riVe 2 



m 



a 2 + /3 2 ) sm9 ((f)i + (f) r 



xE 

it 



72 



(29) 



(30) 



(31) 



uf ko - uj 2 - 2z7 u; ui kQ - uj 2 - 2i7oCJ 
We know that for the complex rotatory power we have the relation p(uj 



(32) 



(ri[ — n r ), further 0; + (fi r = ud (n; + n r ) /c and nf — n 2 = (rii + n r ) (n; — n r ). 
Using these relation we get for p (u) the formula 



2c 



2nNde 2 



-uj 



(a 2 + f3 2 ) sm6 



*E 

k 



-f, 



u 2 k0 - uj 2 - 2i-f u u)l kQ -u 2 - 2ij u 



(33) 
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Now we may find the real and the imaginary part of p(u>) that from ([I]) are 
ORD and CD 



Rep(to) = p{uj) = Atu 2 ■ 



+ 



/<72 (' 



"I 



— UJ 



(< -^ 2 ) +4^2 (u, 2 2 fe0 -u/) + 4 7o 2 ^ 



(34) 



lmp(u) = a{u) = 2A'y u 3 ■ ^ 



-f, 



+ 



'ii 



(<-^ 2 ) +47o 2 ^ 2 « - - 2 ) 



+ 47^ 2 

(35) 

where A = ^^(a 2 + (3 2 ) sin9. 

We know that the quantum transition is optically active when it has a nonze- 
roth value of the rotational strength. It is well known that the formulae for ORD 
and CD are often expressed as the sum of products of rotational strengths for 



different transitions and frequency dependent factors Our formulae (B3) 

and (|35D for ORD and CD do not contain the rotational strengths explicitly. But 



we have showed in the paper JHj that we can express our formulae in terms of 
rotational strengths in a similar way. 

For a transition from some ground state |0) to the excited state \k) the rota- 
tional strength is given by imaginary part of the scalar products of electric and 

magnetic dipole moments (0|d|fc) and (k\m\Q). The rotational strength is so 



R ok = lm((Q\d\k)(k\m\0)) . 



(36) 



In our case of coupled oscillators the ground and the excited state are split 
into two states |^q) and \i]k), rj — 1,2. Therefore the rotational strengths must 
be considered between corresponding ground and excited states and they are 



R 



not, 



Im (vo\d q Jr]k){Vk\rn qjl \r]o) ; t] = 1,2 



(37) 



and for these rotational strengths we have derived in jl4j the results 

he 2 d{a 2 + (3 2 )f qi sm6 



and 



R 2n 



8mc 

he 2 d{a 2 + (3 2 )f q2 sm6 



(38) 



(39) 



We see that all terms of those rotational strengths that depends on the crystal 
structure are comprehended in the formula for the complex rotatory power (j33|) 
and consequentially in the formulae (|3"4]) and (RH) for ORD and CD. 
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3 Discussion 



In the relations (p4|) and (|35|) we restrict only to one optically active quantum 
mechanical transition from the ground states \rjo) to the first excited states |?7i). 
In this case we have for the transition frequencies the relations ujf w = uj\ = uj 2 + Q 
uj 2 -Q. 

(|34]) and (^) can be rewritten in the following way. The results 
for two-oscillator model of optical activity are often expressed as the sum of two 
terms. The first term contains the difference of oscillators strengths f q2 — f qi and 
the second one the sum f qi + f q2 . If we neglect the members containing Q 2 and 
Q'Jq we can write for ORD and CD the relations 



and uj 2 = uj 2 = 
The formulae 



p(uj) = Auj' 



U) 


[uj 2 - uj 2 ) 


+ 


+ U) 


'H - 


uj 2 ) - 




uj 2 ) 2 


+ 4 7o 2 u; 2 




N) - 


^ 2 ) 2 + 


- 47o^ 2 " 


2 



(40) 



a (uj) 



2A70U/ 



fq2 fl 



+ 



2Q (f qi + f q2 ) (uj 2 - uj 2 ) 



(u 2 - u 2 Y + A^uj 2 [ (w g 



UJ 



1 ^ + A 1 2 uj 2 



(41) 



but these results are formally the same as the results for the classical model of 
the optical activity based on the model of two coupled oscillators M . We assume 
that the terms, which Janku has obtained in the quantum mechanical model 
[[| in addition to the classical model, are only the consequence of used solution 
method. 

The results of the solution of the rotational strengths ( |38l) and (^) enable us 
to express the results for p(uj), p(u) and <j(uj), they are given by formulae (p3|), 
(|34"D and (|35D , in another form. The formula for the complex rotatory power can 
be expressed as 



. . 16vr^ 2 2 _ 
P» = ^— EE 



R 



'/OA- 



Tj=l k T]kO 

and similarly the formulae for p(u) and <j(uj) are 



UJ 



p{u) 



criuj 



IQtiNlu 2 2 ( 



UJ 



'IkO 



-uj 2 )R 



L V0k 



he 



V=l k UJ, 



327riVcj 3 7o ^ 



R 



> n >k 



he 



V=l k UJ, 



VkO 



UJ' 



A^Iuj 2 



(42) 



(43) 



(44) 



We can accept these results as the crystal analogs of the similar results which 
we know from the optical activity of molecules [T5|] or recently of delocalized 
molecular aggregates ||16||. In the future it will be verified if these results are 
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valid also for other models of coupled oscillators in the crystal optical activity or 
by what way they must be modified. 
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